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OO I ABSTRACT 

£N| , We propose a zero-point photometric calibration of the data from the ACS/WFC on board the 

CO 1 Hubble Space Telescope, based on a spectrum of Vega and the most up to date in-flight trans- 

mission curves of the camera. This calibration is accurate at the level of a few hundredths of a 
magnitude. The main purpose of this effort is to transform the entire set of evolutionary mod- 
els by Pietrinferni et al. (2004) into a simple observational photometric system for ACS/WFC 
data, and make them available to the astronomical community. We provide the zero points for 
• the most used ACS/WFC bands, and give basic recipes for calibrating both the observed data 

Qh and the models. 

We also present the Colour Magnitude Diagram (CMD) from ACS data of 5 Galactic 
£^ ■ globular clusters, spanning the metallicity range —2.2 <[Fe/H]< —0.04, and provide fiducial 

| points representing their sequences from several magnitudes below the turnoff to the red gi- 

■ ant branch tip. The observed sequences are compared with the models in the newly defined 

photometric system. 



Key words: techniques: photometric - stars: imaging - Hertzsprung-Russell (HR) diagram. 



1 INTRODUCTION 

After two and a half years of operation, the Advanced Camera for 
Surveys (ACS) on board HST has collected a huge number of ob- 
servations of stellar systems, and produced a growing database of 
colour magnitude diagrams (CMDs) and luminosity functions. In 
order to be able to exploit fully the information contained in these 
data sets, and to extract a number of fundamental parameters we 
are interested in, such as metal content, ages, distances, mass func- 
tions, etc, we need to compare the observations with theoretical 
models. The basic requirement for this comparison is the conver- 
sion of stellar models into the ACS photometric observational do- 
main. This conversion is the purpose of the present paper. We have 
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used the most up to date information on the ACS specifications 
to define a natural, in-flight photometric system and to transform 
the Pietrinferni et al. (2004) models into the same system. While 
it is foreseeable that the ACS photometric calibration will improve 
somewhat as more standard fields are studied and cross compared, 
the in-flight transmission efficiency curves are now reasonably well 
established and we would not expect them to change very much. It 
therefore makes sense to provide a conversion of theoretical stellar 
models into the observational domain so that more quantitative and 
qualitative interpretations can be made of observations. 

It is worth mentioning that most astronomical photometric 
investigations are based on some "standard" photometric system. 
Here, by standard, we mean a photometric system that has been 
widely used for a long time in different observatories (e.g., the 
Johnson-Kron-Cousins, Stromgren, Thuan-Gunn, etc systems, see, 
e.g., Landolt-Bornstein 1982), and possibly defined by a number of 
standard stars well distributed in the sky, whose flux has been care- 
fully measured. The calibration to a standard system is sometimes 
the only way to interpret properly the collected data, e.g. when we 
need to compare our photometric sequences with others collected 
with different instruments or at different epochs, or with standard 
results from theoretical models. Still, it must be clearly stated that 
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when the transmission curves of the equipment used to collect the 
observations are rather different from those of any existing stan- 
dard system, the transformation of the data to a standard system is 
difficult, and can be unreliable, particularly for extreme stars (i.e., 
extreme colours, unusual spectral type, high reddening, etc). 

A good example of a problem that one can encounter in using 
even slightly different filter sets at the same telescope is given in 
Momany et al. (2003). 

Many of the ACS filters differ substantially from the "stan- 
dard" filters, and the high photometric precision possible with ACS 
makes the systematic differences more significant. Thankfully, the 
filter and instrumental response are better characterized for ACS 
than they are for the typical ground-based observatory, so that it is 
possible to avoid standard systems altogether and compare the data 
and the model in the same system: the observational domain. All 
that is required for this is to determine an empirical zeropoint for 
each filter by observing a reference star with a well-characterized 
spectrum in order to characterize accurately the throughput and re- 
sponse of the detector for each filter. The theoretical models can 
then be integrated over the same spectral response and compared 
directly with the data. Such a procedure for comparing models with 
data is far preferable to the typical procedure of converting both 
models and data into a "standard" system. 

Therefore, we calculated a set of empirical zero points that 
will allow direct comparison of observational data and theoretical 
models. We have based our calibration on the reference spectrum of 
Vega, and so we will refer to our photometric system as the "ACS 
Vega-mag system". In our derivation of zero points, we used the 
in-flight sensitivity, measured by Sirianni et al. (2002, and Sirianni 
2003, private communication). 

In this paper we have focused our attention on the Wide Field 
Channel (WFC) of the ACS, but the same procedure can be used to 
measure the photometric zero points for the High Resolution Chan- 
nel (HRC). 

In the next Section, we describe briefly the total transmis- 
sion curves, the procedure followed to obtain the zero points, and 
give the aperture corrections. In Section 3, we briefly describe the 
adopted stellar evolution models, our procedure to transform the 
models into the observational plane, and the method used to take 
into account interstellar absorption. In Section 4 we give the zero 
points and the formulae to convert both the theoretical tracks and 
the observed data to the ACS Vega-mag photometric system. Sec- 
tion 5 describes the HST archive data used to test the models, and 
shows, as an example, how these data have been photometrically 
calibrated and compared with our models in the ACS Vega-mag 
plane. Conclusions follow in Section 6 together with instructions 
on how to download the transformed models. 



2 DEFINING THE ACS VEGA-MAG PHOTOMETRIC 
SYSTEM 

In this section, we define our photometric system and will refer to 
it as the ACS(AVFC) Vega-mag system. 

2.1 Total transmission of the HST+ WFC/ACS 

The first step toward the definition of our photometric system is 
to calculate the total transmission curves, T(A) for each individual 
band. These are the products of the individual transmission curves, 
t,(A), for 5 different terms. 

T{X) =fi(A) x ... xt 5 {X) 




4 000 6000 8000 10" 

A[A] 



Figure 1. ACS/WFC transmission curves for the window of the dewar, of 
the WFC/ACS's Mirrors, of the CCD, of the Optical Telescope Assembly 
(OTA), and finally the transmission curves for the two ACS/WFC filters: 
F606W and F814W. The resulting total transmission curve for the two filters 
is shown as a thick line. 

In detail, the t,(A) are: t\ the transmission curve of the window of 
the dewar, ti the transmission curve of the filters, the reflectivity 
of the WFC/ACS mirrors (3 mirrors), t\ an average of the two CCD 
quantum efficiencies 1 and t$ the transmission curve of the Optical 
Telescope Assembly (OTA) of HST 2 . In order to outline our proce- 
dure, we consider two commonly used filters, F606W and F814W. 
We will extend the results to other filters in Section 4. All the in- 
dividual transmission curves are shown in Fig. 1, where two thick 
lines show the total system transmissions for the two considered 
filters. 

The transmission curves measured in the laboratory differ sig- 
nificantly from those found on-orbit (Sirianni et al. 2002). In the 
following we will use only the on-orbit adjusted total transmission 
curves obtained by Sirianni et al. (private communication, paper in 
preparation 2004). The uncertainities for these total transmission 
curves are within 0.5% (De Marchi et al. 2004 3 , Sirianni et al. in 
preparation 2004). 

2.2 The Determination of Zero Points 

We describe below the procedure used to define the zero points of 
the system. 

We have used the spectrum of Vega from ftp:// ftp.stsci.ed u / 

1 These curves are available at the STSCI web-site http:// 
www.stsci.edu /hst /acs /analysis /reference_files/ 
wfc_synphottable_list.html. 

2 This last component has been obtained from the web-site ftp:// 
ftp.stsci.edu /cdbs /cdbs2 /comp /ota/ tabulated in the file 
hst_ota_007_syn.fits. 

3 ISR ACS 2004-08, published after the submission of the present work. 
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cdbs/ cdbs2/ grid/ k93models/ standards/ vega_reference.fits. 
Details of this spectrum can be found at http://www.stsci.edu/ 
instruments/ observatory/ PDF/ scs8.rev.pdf; hereafter it will 



be denoted as F°(A), in units erg s cm _2 A _1 . 

The integrated flux of the assumed spectrum of Vega under the 
total transmission curve T of the filter F606W + instrument is: 



7°(F606W)= I" f F°(X)T } 



'F606W 



(X)dX : 



= 2.7591 x 1(T 6 erg s^'cm" 2 . 

This number is the quantity that theoreticians will need to rescale 
their models to the photometric system we adopt (see Sec. 4). 
The estimated number of expected photo-electrons (e~) is: 

*°(A), 



N°- (F606W) 



A™* F° 



Kim e rW 



r F606w(A)dA ■■ 



= 795,861 e~ r'crn" 2 , 

where the photon energy ey{X) is given by 

ey(A) = hc/X erg. 

Therefore, for Vega, the counts detected by the system in dig- 
ital numbers (DN) would be: 

/g N (F606W) = N°- (F606W) x A eff /GAIN = 

= 3.6004 x 10 ll) /GAIN DN s" 1 , 

where A e g is the effective collecting area of the telescope (45,239 
cm 2 , with the secondary obscuration already accounted for), and 
GAIN is the conversion factor from e~ to DNs. 

At this point, we can define the zero point of our magnitude 
system in the F606W band by imposing mF606W — f° r the mag- 
nitude of Vega (with the assumed spectrum): 



Zp 



F606W _ 



2.5 log 



/ » N (F606W) 
[DNs" 1 ] 



26.391 



For simplicity, we have defined this zero point for GAIN = 1 

e"/DN. 

Since FLT images provided by STSCI are already in photo- 
electrons (e _ ) rather than in DNs (note that DRZ images are given 
in e~ s ) hereafter, to avoid any confusion, we will not consider 
the GAIN, and we will refer to / e - rather than to /qn- 

In an analogous way we can define the zero point of the 
F814Wband: 

Zp F814W = 25 492) 

and in a similar way for all the other filters. 



2.3 Calibration of the Observations 

The formula that observers should use to calibrate their obser- 
vations into our WFC/ACS Vega-mag system magnitudes in the 
F606W band is: 

" -^Aperture 



«F606W = -2.5 log 10 



exptime 



-Zp 



F606W 



where exptime is the exposure time, and ZAperture is the aperture 
correction. In the next section we discuss in detail the aperture 
correction, which we will always treat as a positive quantity to be 



JTTTTJ7P 


A>»AP(~)-AP(~) 


^ m AP(0".5)-AP(=o) 


Am AP(0".3)- 


F435W 


0.000 


0.095 


0.132 


F475W 


0.000 


0.084 


0.120 


F502N 


0.000 


0.078 


0.114 


F550M 


0.000 


0.076 


0.110 


F555W 


0.000 


0.076 


0.110 


F606W 


0.000 


0.079 


0.115 


F625W 


0.000 


0.080 


0.118 


F658N 


0.000 


0.079 


0.123 


F660N 


0.000 


0.079 


0.123 


F775W 


0.000 


0.078 


0.123 


F814W 


0.000 


0.079 


0.130 


Table 1. Am 


AP(r)-AP(-o) corrections to be added to the Zero Points. 



AP(°< 



subtracted. This calibration is intended for the officially calibrated 
DRZ images which are already corrected for: flat field, dark, bias, 
and geometrical distortion. 4 



2.4 Aperture Corrections 

Breathing of the telescope tube, the complex optical system, and 
charge diffusion induce temporal and spatial variation of the core 
of the point-spread function (PSF), such that the 1 st -ring ratio of 
the diffraction figure changes considerably across the chip and at a 
minor level with orbital phase. 

For these reasons, we choose to refer our zero points to an 
infinite aperture. That makes it easier for those studying surface 
photometry of extended sources to calibrate their magnitudes to 
this system. The aperture corrections for extended sources are null, 
being the Zero Points of the system referred to an infinite aperture: 5 



perture 



-0. 



On the other hand, in practice all the energy of a point source is 
contained within 4" (~ 80 pixels), and ~ 93% within 0."5 (~ 10 
pixels). In dense environments, where crowding makes an aperture 
correction to an infinite aperture impossible, one has to rely on the 
in-flight encircled energy. 6 This means that one has to correct the 
instrumental magnitude (for example the PSF-fitted magnitude) to 
the magnitude at a specified aperture r by adding a Am PSF _ AP ( r ), 
and then sum an additional quantity AmAp( r )_AP(<>o)> obtained from 
the tabulated encircled energy in order to link the magnitudes to the 
infinite aperture for which we have calculated the Zero Points. 
In summary, 



^Aperture : 



- Am, 



PSF-AP(r) 



+ Am 



AP(r)-AP(» 



In Table 1 we give the AmAp( r )_AP(oo) corrections to be added 
to the zero points for two circular apertures of radius 0."3 (~6 pix- 
els) and 0."5 (as in Holtzman et al. 1995) and for the different fil- 
ters. For clarity, we also report in Tab. 1 the correction to zero points 
for an infinite aperture, which is zero by definition. 

For point sources, the Amp SF _ A p( r ) correction should be esti- 
mated in each individual image, remembering that the smaller the 



4 The images are still affected, at least by CTE, and charge diffusion (see 
Riess & Mack, ISR 2004-06, and Krist, ISR 2003-06). In Section 5. 1 we use 
FLT images, and describe how we linked our photometry for those images 
to the officially corrected DRZ images. 

5 Note, that this is really only true for infinite extended sources, or sources 
larger than the widest part of the PSF halo. 

6 This is tabulated at the STSCI web-site: http:// www.stsci.edu 
/hst/acs/analysis/ reference_files/ wfc_synphottableJist.html 
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aperture (r), the more the photometry is affected by the problems 
mentioned at the beginning of this section (i.e. spatial variability of 
the PSF and charge diffusion). See Section 5.2 for uncertainties on 
the values of Am PSF _ A p(,.) . 

Moreover, there appears to be a second-order effect which de- 
pends on the colour of the stars and on the filter used (Sirianni, pri- 
vate communication, paper in preparation). In practice, the present 
calibration is useful only if the degree of accuracy desired by the 
observers is no better than 2-3 hundredths of a magnitude. How- 
ever, taking into account all the problems of ACS, it seems hard to 
achieve a much higher accuracy. It should be noted that the pho- 
tometric stability of WFPC2 was certainly not any better than this, 
but the higher precision of WFC/ACS makes systematic errors (on 
the zero points) the dominant source of error in the final photome- 
try. 

We have decided not to deal with the two filters: F850LP and 
F892N, because they are much more problematic (see Sirianni, in 
preparation). 



3 TRANSFORMING THE STELLAR MODELS INTO 
THE ACS VEGA-MAG SYSTEM 

3.1 The Stellar Evolution Models 

In this work we adopt the stellar-evolution model database in 
Pietrinferni et al. (2004); for a detailed description of the models, 
as well as a comparison with other stellar model libraries available 
in the literature, we refer the interested reader to this reference. For 
the purposes of this paper, we briefly list the main characteristics 
of this model library. 

Evolutionary tracks are available for up to 41 mass values at 
each of 10 selected initial chemical compositions. The minimum 
mass considered is roughly 0.5M©, while the maximum value is 
always equal to IOMq. The whole set of evolutionary computa- 
tions has been used to compute isochrones for a large age range, 
from ~ 40 Myr to ~ 14 Gyr, covering all stellar evolution phases 
from the Zero Age Main Sequence up to the first thermal pulse on 
the Asymptotic Giant Branch or to C ignition. For each chemical 
composition, a Zero Age Horizontal Branch (ZAHB) locus, as well 
as post-ZAHB evolutionary tracks were also calculated for a large 
range of masses, by employing the He-core mass and chemical pro- 
files in the H-rich envelope of the Red Giant Branch (RGB) progen- 
itor having an age at the He-flash of the order of ~ 12—13 Gyr. 

The set of models has been computed for metallicities in the 
range 0.0001 ^ Z ^ 0.04 for both a scaled-solar mixture and an 
a-enhanced mixture (<[a/Fe]>=0.4; paper in preparation). 

The adopted reference scaled-solar heavy element mixture is 
from Grevesse & Noels (1993). As for the He-abundance, the mod- 
els employ a primordial mass fraction of lo=0.245 , based on the 
estimates by Cassisi et al. (2003, see also Salaris et al. 2004) ob- 
tained from the ^-parameter method applied to a large sample of 
Galactic globular clusters (GGCs). This Yq value is in fair agree- 
ment with recent independent determinations of the cosmological 
baryon density from the WMAP results (Spergel et al. 2003). To re- 
produce the initial solar He-abundance obtained from an appropri- 
ately calibrated solar model, an He enrichment ratio dY jdZ = 1 .4 
has been used (Pietrinferni et al. 2004). 

All models have been computed with outer boundary con- 
ditions obtained by integrating the atmospheric layers with the 
Krishna-Swamy (1966) T(t) relationship. Superadiabatic convec- 
tion is treated according to the Cox & Giuli (1968) formalism of 



mixing-length theory (Bohm-Vitense 1958), and the mixing-length 
parameter has been fixed by the solar calibration (ml = 1.713) and 
kept constant for all masses during all evolutionary phases; mass 
loss is included by employing the Reimers formula (Reimers 1975) 
with the free parameter 77 set to 0.4. 

At each chemical composition, models have been computed 
with and without overshooting from the convective cores during the 
main sequence (MS) phase. In the former case, Xov = 0.20 x Hp for 
masses larger than 1.7M , where Xov is the extension of the over- 
shooting region in units of the local pressure scale height. For stars 
less massive than 1.1M Q , Xov = 0, while in the intermediate range 
Xov varies according to the relation Xov~(M/M e -0.9)/4. This 
choice allows a smooth variation of the isochrone turn-off (TO) 
morphology, and a smooth decrease to zero of the convective core 
mass. This parametrization also provides a good match to the TO 
morphology in the CMD of Galactic open clusters of different ages. 
Induced overshooting and semi-convection during the central He- 
burning phase are accounted for following Castellani et al. (1985). 

Radiative opacities from the OPAL tables (Iglesias & 
Rogers 1996) for temperatures larger than 10 4 K, and from Alexan- 
der & Ferguson (1994) for lower temperatures have been em- 
ployed, supplemented by opacities for electron degenerate mat- 
ter computed following Potekhin (1999). The relevant energy loss 
rate from plasma-neutrino processes comes from Haft, Raffelt & 
Weiss (1994), whereas for all other processes the reader is referred 
to Cassisi & Salaris (1997); the nuclear reaction rates are from 
the NACRE database (Angulo et al. 1999), with the exception of 
the 12 C(a,y) 16 reaction, in which case the more accurate recent 
determination by Kunz et al. (2002) has been adopted. Electron 
screening is treated according to Graboske et al. (1973). 

The accurate Equation of State (EOS) by A. Irwin has been 
used. An exhaustive description of this EOS is in preparation (Ir- 
win et al. 2004) but a brief discussion of its main characteristics can 
be found in Cassisi et al. (2003). This EOS covers the full stellar 
structure for the whole evolutionary phase and mass range spanned 
by the model database; in Pietrinferni et al. (2004) and Bahcall, 
Serenelli & Pinsonneault (2004) it is shown that this EOS and the 
latest version of the OPAL EOS (Rogers 2001) provide almost iden- 
tical results in the case of the Sun. The accuracy and reliability of 
the present evolutionary scenario has been successfully tested in 
Cassisi et al. (2003), Riello et al. (2003) and Salaris et al. (2004) 
in the regime of Galactic globular clusters, and by Pietrinferni et 
al. (2004) in the regime of field stars and open clusters. Cassisi 
et al. (2003) and Pietrinferni et al. (2004) have also compared our 
adopted set of theoretical isochrones with other existing isochrone 
databases widely cited in the current literature. This comparison 
provides an estimate of the present uncertainties in stellar evolu- 
tion models due to different choices for the stellar input physics 
(see also Chaboyer 1995). The main conclusions of this exercise 
are summarized below. 

Isochrones for typical globular cluster ages (of the order of 
10 Gyr) show a remarkable consistency for the prediction of the 
Main Sequence location and Turn Off luminosity at a given ini- 
tial chemical composition, whereas the Red Giant Branch effec- 
tive temperature and Zero Age Horizontal Branch brightness dis- 
play values that differ by amounts of ~ 100-200 K and ~0. 1 mag, 
respectively. Globular cluster ages estimated from the difference 
between Turn Off and Horizontal Branch brightness would con- 
sequently be uncertain by ~10 % (Chaboyer 1995 determined a 
similar uncertainty, of the order of 15%). 

Isochrones for ages lower than ~4-5 Gyr show again con- 
sistent Main Sequence locations, but differing Turn Off luminosi- 
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ties, effective temperatures, and Helium burning luminosities, due 
mainly to the uncertain treatment of the overshooting from the 
Main Sequence stars' convective cores. These uncertainties cause 
variations in the estimated age of stellar populations within this age 
range, of at most ~30 %, Also in this age range red giants' effective 
temperatures show a spread of ~ 100-200 K. 



3.2 The Transformation to the Observational Plane 

Here we describe the procedure we have followed to transform our 
models into the ACS Vega-mag system defined in Sect. 2. 

We computed ACS colour indices by using the homogeneous 
set of ODFNEW model atmospheres and synthetic fluxes 7 (Castelli 
& Kurucz 2003) computed with the Kurucz ATLAS9 code. We gen- 
erated grids of indices for the following values of [Fe/H]: 0.5, 0.2, 
0.0, -0.5, -1.0, -1.5, -2.0, -2.5. 

In each grid point, T e ff ranges from 3500 to 50,000 K, logg 
from 0.0 to 5.0 and the microturbulent velocity is ^=2.0kms _l . 
We tabulated for each T e g and logg gridpoint the BCy bolometric 
correction, the V magnitude, and the colour indices V — mF435W> 
V - mF475W- V - m F5 55 W , V - m F606w , V - m F 625W, V - m F 7 75w , 
and V — >«F814W- The ACS magnitudes were computed by using 
the WFC/ACS transmission curves described in Sect. 2.1, while 
we adopted the V passband from Bessel (1990). 

We adopted V = 0.03 for Vega and set all colour indices (V — 
m filte?) to zero, so that the Vega ACS magnitudes would be equal 
to 0.00 in all passbands. The Vega model and flux are the ATLAS9 
model and the ATLAS9 flux from Castelli & Kurucz (1994) with 
parameters r eff =9550K, logg=3.95, [M/H] = -0.5, §=2kms _1 . 
Note that this model has been built specifically to reproduce the 
observations of the spectrum used in Sect. 2.2. The absolute zero- 
point for the V magnitude was fixed by adding the constant —21.10 
mag to the computed V for Vega (Bessell, Castelli & Plez 1998). 

In Table 2 we show an example of the theoretical isochrones 
transformed into the ACS Vega-mag photometric system. The 
header of the isochrone table provides information about the num- 
ber of lines in the file (Np), the total metallicity [M/H], the abun- 
dance by mass of metals (Z) and helium (Y), and the correspond- 
ing age in billions of years, t(Gyr). The content of the columns 
in the following lines is the following (from left to right): 1) the 
initial mass of each star, 2) the current mass at the specified age 
(computed accounting for mass loss - see previous section), 3) the 
logarithm of the surface luminosity in solar units, 4) the logarithm 
of the effective temperature, 5) the ACS MF435W magnitude, 6) the 
ACS Mp475w magnitude, 7) the ACS A/F555W magnitude, 8) the 
ACS Mp 6( )6 W magnitude, 9) the ACS M F g25w magnitude, 10) the 
ACS Mp775w magnitude, 11) the ACS Mpgi4w magnitude. 

The entire set of transformed isochrones, as well as the ZAHB 
and post-ZAHB evolutionary tracks can be retrieved from the fol- 
lowing URL: http:// www.te.astro.it/BASTI/index.php. 



fluxes and the reddened WFC/ACS magnitudes. The Agi ter extinc- 
tions in the WFC/ACS filters were obtained as differences between 
the reddened and unreddened m^iter magnitudes. 

Tables 3 and 4 give the values of the Agi ter for a cool star and 
a hot star respectively. We give these values for the main filters, 
F435W, F475W, F555W, F606W, F625W, F775W, and F814W. We 
note how these values are very close to those tabulated by Holtzman 
et al. (1995) for WFPC2's similar filters. 



4 ACS VEGA-MAG SYSTEM ZERO POINTS 

Table 5 lists the zero points for the ACS Vega-mag system in other 
useful pass-bands; the zero points for all the other filters can be 
easily obtained following the recipes outlined above. 

We note that this calibration is linked to the chosen spectrum 
of Vega. These zero points are probably accurate to 2-3 hundredths 
of a magnitude, and can (and hopefully will soon) be updated, once 
more accurate data on the ACS performance is available. 

The main sources of systematic error in the calibration of the 
observations are: residuals in correction for CTE, flat fielding, aper- 
ture corrections, charge diffusion, and the encircled energy which 
seems to depend on stellar colour (Sirianni, private communica- 
tion). Still, the lack of any other official and complete photometric 
calibration of the ACS, useful for the comparison with stellar mod- 
els, has led us to calculate and publish this calibration. 

In summary, photometric observations can be transformed into 
the ACS/WFC Vega-mag system defined in this paper by applying 
the following relation: 



^filter = -2.5 log 



L0 



exptime 



-Zp filter + 



— Am', 



filter 

PSF-AP(r; 



— Am' 



filter 

AP(r)-AP(o 



(1) 



Where r can assume the values: °°, 0".5, or 0".3. Note that the last 
two positive quantities in Eq.Qare subtracted. 

Models can be transformed into the observational plane of the 
same system by applying: 



^filter = -2.5 log 
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/(filter) 
/"(filter) ' 



(2) 



Here the meaning of /(filter), and /"(filter) are clear from the dis- 
cussions in previous sections. In Sect. 5.2 we give a practical ex- 
ample of how to calibrate the observations. 

Our Zero points agree well with those presented in the recent 
ISR by De Marchi et al. (2004-08). 8 The only differences are due 
to a different choice of spectrum of Vega (Sirianni private comuni- 
cation). Nevertheless, for broad filters the Zero points agree within 
~0.005 mag. 



3.3 Reddening in the WFC/ACS Vega-mag system 

Following the recipe and the motivations given by Holtzman et 
al. (1995) we also computed the interstellar absorption in the 
WFC/ACS bands as a function of E(B — V). By using the extinction 
curve from Table 1 in Mathis (1990) for Ry — 3.1 we computed 
for different values of E(B — V) the interstellar reddening at 
the wavelengths of the ATLAS9 flux, the corresponding reddened 



5 COMPARISON OF MODELS WITH OBSERVATIONS 
5.1 Observations and Data Reduction 

In order to give a practical example of our procedures and to pro- 
vide a useful test for our models, we used archival WFC observa- 
tions of 5 Galactic globular clusters (GO-9453, PI: Brown), which 
span a wide range in metallicity (—2.2 < [Fe/H] < —0.04). The 



7 Available at http:// wwwuser.oat.ts.astro.it/ castelli/ grids.html. 



1 http:// www.stsci.edu/ hst/ acs/ documents/ isrs/ isr0408.pdf. 
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Isochrone by Pietrinfemi - Cassisi - Salaris - Castelli 2004 Standard Model 



Scaled solar model & transformations for ACS (Castelli 2004) 
Afp = 2000, [M/H] = —0.659, Z = 0.0040, 7 = 0.251, t(Gyr) = 11.0000 



(M/M e ) in 


(M/Mq ) 


logL/L 


logT cff 


M F435W 


M F475W 


M F555W 


M F606W 


M F625W 


M F775W 


M F814W 


.5000000000 


.4999030025 


-1.20676 


3.62790 


9.540 


9.063 


8.418 


8.024 


7.693 


7.069 


6.982 


.5011654487 


.501 1233468 


-1.20259 


3.62829 


9.525 


9.048 


8.405 


8.012 


7.681 


7.059 


6.972 


.5023308974 


.5023308974 


-1.19841 


3.62868 


9.510 


9.034 


8.392 


8.000 


7.670 


7.049 


6.962 


.9066248329 


.6047402926 


2.89310 


3.60793 


-0.175 


-0.871 


-1.635 


-2.045 


-2.382 


-3.053 


-3.153 


.9066266265 


.6037663031 


2.92275 


3.60503 


-0.190 


-0.899 


-1.675 


-2.093 


-2.435 


-3.116 


-3.218 


.9066284200 


.6026715201 


2.95336 


3.60202 


-0.204 


-0.926 


-1.715 


-2.143 


-2.490 


-3.181 


-3.285 



Table 2. Example of an isochrone available at http:// www.te.astro.it/ BASTI/ index. php. 



T sff = 4.000K, logs = 4.50, [M] = 0.00, v turb = 2.00, l/H = 1.25 
E(B-V) ^F435W - 4 F475W A F555W A F606W A F625W ^F775W A F814W 



0.00 0.000 

0.05 0.202 

0.10 0.404 

0.15 0.606 

0.20 0.808 

0.25 1.010 

0.30 1.212 

0.35 1.413 

0.40 1.614 

0.45 1.815 

0.50 2.017 

0.55 2.217 

0.60 2.418 

0.65 2.619 

0.70 2.819 

0.75 3.020 

0.80 3.220 

1.00 4.020 

1.25 5.016 

1.50 6.010 

1.75 7.001 

2.00 7.990 

2.25 8.975 

2.50 9.959 

2.75 10.940 

3.00 11.919 

3.25 12.896 

3.50 13.871 

3.75 14.844 

4.00 15.815 

4.25 16.785 

4.50 17.753 

4.75 18.719 

5.00 19.684 



0.000 0.000 

0.182 0.158 

0.365 0.316 

0.546 0.474 

0.728 0.632 

0.909 0.790 

1 .090 0.947 

1.271 1.104 

1.451 1.262 

1.631 1.419 

1.811 1.575 

I. 990 1.732 
2.169 1.889 
2.348 2.045 
2.527 2.201 
2.705 2.357 
2.883 2.513 
3.592 3.136 
4.472 3.911 
5.347 4.682 
6.215 5.450 
7.077 6.216 
7.935 6.978 
8.788 7.737 
9.636 8.494 
10.480 9.249 

II. 320 10.001 
12.157 10.751 
12.991 11.499 
13.822 12.244 
14.650 12.988 
15.476 13.730 
16.299 14.470 
17.120 15.209 



0.000 0.000 

0.140 0.132 

0.279 0.263 

0.419 0.394 

0.558 0.526 

0.696 0.657 

0.834 0.788 

0.972 0.919 

1.110 1.050 

1.248 1.180 

1.385 1.311 

1.522 1.442 

1.658 1.572 

1.794 1.702 

I. 930 1.833 
2.066 1.963 
2.202 2.093 
2.741 2.612 
3.409 3.260 
4.072 3.905 
4.730 4.547 
5.382 5.188 
6.031 5.827 
6.675 6.463 
7.315 7.098 
7.951 7.731 
8.585 8.362 
9.215 8.991 
9.843 9.619 
10.468 10.245 

II. 091 10.870 
11.711 11.493 
12.329 12.115 
12.946 12.735 



0.000 0.000 

0.099 0.092 

0.198 0.184 

0.297 0.276 

0.396 0.367 

0.495 0.459 

0.594 0.550 

0.692 0.641 

0.791 0.732 

0.890 0.823 

0.988 0.913 

1 .086 1 .004 

1.185 1.094 

1.283 1.184 

1.381 1.274 

1.479 1.364 

1.577 1.454 

1.968 1.811 

2.454 2.254 

2.939 2.693 

3.422 3.128 

3.902 3.559 

4.381 3.987 

4.858 4.411 

5.333 4.832 

5.806 5.250 

6.277 5.664 

6.747 6.075 

7.215 6.483 

7.681 6.888 

8.145 7.291 

8.608 7.690 

9.070 8.088 

9.530 8.483 



Table 3. Extinctions (in mag) in different ACS filters as a function of the colour excess E(B-V) of a cool star. 



r efr = 40, 0O0K, logg = 4.50, [M] = 0.00, v turb = 2.00, l/H = 1.25 



E(B-V) 


/1 F435W 


/1 F475W 


A F555W 


"F606W 


A F625W 


J F775W 


^FSMW 


0.00 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.05 


0.212 


0.194 


0.165 


0.153 


0.135 


0.102 


0.096 


0.10 


0.423 


0.387 


0.329 


0.305 


0.269 


0.203 


0.192 


0.15 


0.635 


0.580 


0.493 


0.456 


0.403 


0.304 


0.288 


0.20 


0.846 


0.773 


0.657 


0.608 


0.538 


0.406 


0.383 


0.25 


1.057 


0.965 


0.821 


0.759 


0.672 


0.507 


0.479 


0.30 


1.267 


1.157 


0.985 


0.909 


0.806 


0.608 


0.574 


0.35 


1.478 


1.348 


1.148 


1.059 


0.940 


0.709 


0.669 


0.40 


1.688 


1.539 


1.311 


1.209 


1.073 


0.810 


0.764 


0.45 


1.898 


1.730 


1.474 


1.358 


1.207 


0.911 


0.859 


0.50 


2.108 


1.921 


1.637 


1.507 


1.341 


1.012 


0.953 


0.55 


2.318 


2.111 


1.800 


1.656 


1.474 


1.112 


1.048 


0.60 


2.527 


2.300 


1.963 


1.804 


1.608 


1.213 


1.142 


0.65 


2.737 


2.490 


2.125 


1.951 


1.741 


1.313 


1.237 


0.70 


2.946 


2.678 


2.288 


2.099 


1.874 


1.414 


1.331 


0.75 


3.155 


2.867 


2.450 


2.246 


2.007 


1.514 


1.424 


0.80 


3.364 


3.055 


2.612 


2.392 


2.140 


1.615 


1.518 


1.00 


4.196 


3.804 


3.258 


2.975 


2.671 


2.015 


1.891 


1.25 


5.233 


4.732 


4.062 


3.694 


3.333 


2.513 


2.354 


1.50 


6.265 


5.651 


4.862 


4.405 


3.992 


3.009 


2.813 


1.75 


7.292 


6.562 


5.658 


5.107 


4.648 


3.503 


3.268 


2.00 


8.315 


7.465 


6.450 


5.800 


5.302 


3.995 


3.719 


2.25 


9.333 


8.360 


7.239 


6.486 


5.954 


4.485 


4.166 


2.50 


10.347 


9.247 


8.023 


7.165 


6.603 


4.972 


4.609 


2.75 


11.357 


10.128 


8.804 


7.837 


7.250 


5.458 


5.048 


3.00 


12.362 


11.003 


9.581 


8.503 


7.895 


5.942 


5.483 


3.25 


13.364 


11.872 


10.355 


9.163 


8.538 


6.423 


5.915 


3.50 


14.362 


12.735 


11.126 


9.818 


9.178 


6.902 


6.343 


3.75 


15.357 


13.594 


11.894 


10.468 


9.817 


7.380 


6.767 


4.00 


16.348 


14.447 


12.658 


11.114 


10.453 


7.856 


7.187 


4.25 


17.336 


15.297 


13.420 


11.755 


11.088 


8.329 


7.604 


4.50 


18.321 


16.142 


14.179 


12.393 


11.721 


8.801 


8.018 


4.75 


19.303 


16.985 


14.935 


13.028 


12.352 


9.271 


8.429 


5.00 


20.283 


17.823 


15.689 


13.659 


12.981 


9.740 


8.837 



Table 4. As in Table 3, but for a hot star. 
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filter 


F435W 


F475W 


F502N 


F550M 


F555W 


F606W 


Zp filter 

7° (filter) 


25.785 
2.0994 


26.168 
2.7611 


22.338 
0.075320 


24.861 
0.69396 


25.718 
1.6028 


26.391 
2.7591 


filter 


F625W 


F658N 


F660N 


F775W 


F814W 




Zp mta 


25.722 


22.381 


21.342 


25.254 


25.492 





Afilter) 1.3687 0.059770 0.022895 0.72740 0.87198 
Table 5. Natural ACS-VEGA System Zero Points (Zp flltel ), and Fluxes 7° (filter) (in 10~ 6 erg s -1 cm" 2 ). 



data are summarized in Table 6. They consist of 1 short, 1 interme- 
diate, and 1 (relatively) long exposure, in both F606W and F814W 
filters, for each cluster. 

We carried out the photometry with algorithms based on the 
effective point-spread function (ePSF) fitting procedure described 
by Anderson & King (2000). We obtained the ePSF from a set of 
images from a different project using the same filters, since the 
dithering of GO-9453 images was unsuitable to constrain the ePSF 
shape (Anderson & King, 2000). 

The resampling present in the DRZ images provided by 
STSCI does not allow the use of the ePSF method, limiting the 
precision of the output photometry. We therefore performed our 
fitting photometry using FLT images. These images, affected by 
geometrical distortion, have been flat-fielded to preserve surface 
brightness, but not flux. Therefore, we must multiply each pixel by 
its projected area to obtain a proper flux for each star. 

We find slight differences of flux-scaling between our cor- 
rected FLT images and the corresponding DRZ images (on the 
order of 1% or so). Since it is the DRZ images that have been 
carefully calibrated via the STSCI pipe-line, we should convert 
any photometry measured on FLT images into the corresponding 
infinite-aperture fluxes in the DRZ images. The photometry has 
been calibrated as described in the previous sections. 

5.2 An Example: Calibration of the 47Tuc Photometry 

For the GO-9453 47Tuc data, the zero point of the long expo- 
sures has been used to calibrate the intermediate and short ones 
in both F606W (for which the exposure time is 70s) and F814W 
filters (for which the exposure time is 72s, see Table 3). We de- 
termined the zero-point difference between our PSF-fitting pho- 
tometry and aperture photometry with a 0".5 aperture for each 
filter. The differences are Am™ 6 *^,, „ = +0.15(±0.015), and 

Am PSF^AP(0" 5) = +0.25(±0.015). To these we need to add the 
aperture corrections. According to Table 1: Am^^^^ _ A p^ = 



Cluster 



F606W data 



F814Wdata 



Am 



F814W 



AP(0".5)-AP(» 



+0.079 (equal for the two filters, at this aper- 



tures). In summary, from Eq. Q: 

70 



'"F606W = -2.5 log 10 



n?F8l4W = -2.5 log 10 



72 



26.391-0.15-0.079 



-25.492-0.25-0.079 



[Note: 1) —2.5 log 1() / e - is simply the instrumental magnitude 
(some codes, such as DAOPHOT, add a constant value of 25 to the 
instrumental magnitude to avoid negative numbers); 2) the adopted 
normalization of the fitting magnitude is taken into account by the 
aperture correction Am PSF _ AP ( r j; 3) the sign of the last two quanti- 
ties in each equation is negative.] 

The calibrated CMDs are presented in Fig. 2. Superposed on 



NGC 6528 
NGC 5927 
NGC 104 
NGC 6752 
NGC 6341 



4s, 50s, 450s 
2s, 30s, 500s 
0.5s, 6s, 70s 
0.5s, 4s, 40s 
0.5s, 5s, 90s 



Is, 20s, 350s 
0.7s, 15s, 340s 
0.5s, 5.5s, 72s 

0.5s, 4s, 45s 
0.5s, 6s, 100s 



Table 6. Data set of the observed clusters. 



the observed CMD are the fiducial lines of the MS, SGB, and RGB. 
The fiducial points which define them have been determined by eye 
and are given in Table 7. Note that, due to the differential redden- 
ing in the field, the fiducial lines in the cases of NGC 6528 and 
NGC 5927 are poorly defined. 

These fiducial points can be used as a reference to com- 
pare against those measured for other globular clusters or galaxies 
whose stars are resolved (Brown et al. 2003). 



5.3 Comparison with the Models 

In this section we compare the theoretical models transformed to 
the ACS Vega-mag system with observed CMDs. 

Figure 3 displays our adopted a-enhanced isochrones and 
ZAHBs compared with the observed CMDs of the GGCs 47 Tuc 
(NGC 104), NGC 6752, and M92 (NGC 6341). These clusters 
span a wide metallicity range (—2.2 < [Fe/H] < —0.6). We selected 
these clusters for the comparison because of their very low extinc- 
tion and lack of differential reddening. 

In the following discussion, in order to make the compari- 
son with the values in the literature more straightforward, we give 
the reddening E(B — V) and the apparent distance modulus in the 
V band (yy). The E[B — V) comes from E(mpa)6W — m F8l4w) = 
^F606W — ^F814W as adopted in the best fitting by interpolation of 
the values in Table 3. The yy comes from the best fitting distance 
modulus in the F606W band (y F 606w): yv =yF606W -^F606W +Av, 
where A v = 3 AE{B-V). 

In the case of 47 Tuc we used a-enhanced (< [a/Fe] >= 0.4) 
isochrones with [Fe/H] = —0.66 (Z=0.008), a value in very good 
agreement with spectroscopic measurements (see, e.g., Percival et 
al. 2002, Gratton et al. 2003); the best match is obtained for an 
apparent distance modulus yy = 13.45, an age of 10 Gyr, and a 
reddening E(B — V) =0.04 (the 8 and 13 Gyr isochrones are also 
shown). These values of distance modulus and reddening are con- 
sistent with empirical determinations by Percival et al. (2002) and 
Gratton et al. (2003). 

The quality of the fit is satisfactory along most of the 
isochrone, worsening at the redder end of both the RGB and MS. 
The deviation of the lower part of the sequence from the models 
for intermediate-metallicity clusters is well known (see Bedin et al. 
2001). 
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For NGC 6752 we have employed our isochrones with 
[Fe/H] = -1.26 (Z=0.002), a value close to the estimate by Rut- 
ledge et al. (1997) with their calibration of the Call triplet index on 
the GGC Carretta & Gratton (1997) metallicity scale. We obtain a 
best match with yy=13.25, an age of 12 Gyr and E(B - V)=0.04 
(the 10 and 14 Gyr isochrones are also shown). 

Again, the reddening agrees with independent empirical es- 
timates (e.g. Gratton et al. 2003), whereas the distance is slightly 
longer than the MS-fitting estimates (Carretta et al. 2000, Gratton 
et al. 2003). These empirical estimates employ a metallicity ~0.2 
dex lower than our adopted value; it is interesting to note, judging 
from the numerical experiment performed by Carretta et al. (2000, 
their Fig. 5), that their empirical MS-fitting distance would come 
into agreement with our value if our assumed metallicity is used in 
their MS-fitting technique. 

Had we used our isochrones for [Fe/H] = —1.57 (Z=0.001) 



and a value close to the metallicity estimate by Rutledge et al. 
(1997) obtained calibrating the calcium index on the Zinn & West 
(1984) metallicity scale, the quality of the fit to the RGB location 
would be worse (isochrone too blue), the age older by ~1 Gyr, 
E(B — V) ~0.02 mag higher and yy ~0.07 mag larger. 

The comparison with the M 92 (NGC 6341) CMD has been 
performed using isochrones with [Fe/H]~ —2.1 (Z=0.0003), a 
value in fair agreement with the estimate by Carretta & Grat- 
ton (1997) ([Fe/H] = —2.16). The best agreement is achieved for 
an age of 12 Gyr, E(B - V)=0.04 and y v = 14.85 (the 10 and 14 
Gyr isochrones are also shown). 

6 SUMMARY 

In this paper we have proposed a simple zero-point photometric cal- 
ibration of the ACS/WFC data, based on a model spectrum of Vega 
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Colour 




















2.57 


12.17 






NGC 6341 






NGC 5927 
Colour Mug 


2.38 
2.05 
1.92 


11.83 
11.39 
11.30 


NGC 6752 
Colour M;ig 


Colour 
1.08 


Mag 
11.67 






2.70 


14.63 


1.75 


11.23 




10.48 
10.99 
11.32 
11.80 
12.42 
13.19 
13.88 
14.60 
15.10 
15.60 
15.91 
16.15 
16.35 
16.44 
16.51 
16.66 
16.78 
16.87 
17.02 
17.16 
17.33 
17.61 
17.83 
18.26 
18.67 
19.17 
19.93 
20.63 
21.56 
22.61 
23.31 


0.98 


12.06 


NGC 6528 
Colour Mug 


2.59 
2.33 
2.03 


14.45 
14.17 
14.15 


1.56 
1.40 
1.27 


11.25 
11.36 
11.55 


1.22 
1.08 
1.00 
0.94 
0.87 
0.81 
0.77 
0.73 
0.71 
0.69 
0.68 
0.66 
0.65 
0.64 
0.61 
0.58 
0.55 
0.53 
0.52 
0.51 
0.51 
0.52 
0.53 
0.56 
0.60 
0.66 
0.78 
0.93 
1.12 
1.27 
1.33 


0.90 
0.80 
0.74 


12.70 
13.73 
14.77 


3.79 
2.91 
2.33 
1.72 
137 


16.46 
15.52 
15.24 
15.36 
15.75 
16.26 
16.85 
17.59 
18.41 
19.33 
19.82 
19.97 
19.97 
19.89 
20.02 
20.35 
20.86 
21.58 
22.45 
23.40 
24.14 
24.88 
25.88 


1.87 


14.18 


1.16 


11.78 


0.72 


15.41 


1.72 
1.57 


14.38 
14.69 


1.07 
1.01 


12.03 
12.35 


0.69 
0.68 


16.06 
16.54 


1.46 


15.02 


0.94 


12.84 


0.67 


16.86 


1.39 


15.40 


0.88 


13.32 


0.66 


17.07 


1.31 


15.93 


0.83 


13.83 


0.65 


17.22 


1.48 
1.43 
1.38 
1.33 
1.31 
1.28 
1.26 
1.20 
1.14 
1.09 


1.24 


16.59 


0.80 


14.26 


0.63 


17.34 


1.18 


17.18 


0.76 


14.90 


0.60 


17.55 


1.15 


17.84 


0.74 


15.55 


0.57 


17.68 


1.12 


18.53 


0.73 


15.98 


0.55 


17.76 


1.12 


18.94 


0.71 


16.37 


0.52 


17.91 


1.11 


19.21 


0.71 


16.58 


0.49 


18.03 


1.10 


19.32 


0.70 


16.74 


0.48 


18.12 


1.08 


19.44 


0.69 


16.81 


0.46 


18.21 


1.04 


19.49 


0.67 


16.90 


0.45 


18.30 


1.01 


19.52 


0.65 


16.95 


0.45 


18.42 


1.08 
1.12 
1.23 
1.36 
1.60 
1.85 
2.10 
2.35 


0.97 


19.64 


0.62 


16.99 


0.45 


18.58 


0.95 


19.75 


0.59 


17.06 


0.45 


18.74 


0.95 


20.00 


0.57 


17.16 


0.46 


18.93 


0.95 


20.33 


0.56 


17.30 


0.48 


19.18 


0.98 


20.66 


0.56 


17.52 


0.49 


19.43 


1.01 


21.09 


0.56 


17.76 


0.52 


19.75 


1.08 


21.63 


0.57 


18.01 


0.55 


20.08 


1.23 


22.44 


0.59 


18.33 


0.60 


20.44 






1.51 


23.58 


0.62 


18.67 


0.66 


20.90 






1.77 


24.47 


0.67 


19.07 


0.72 


21.29 






1.99 


25.36 


0.73 


19.50 


0.81 


21.85 






2.21 


26.18 


0.82 
1.02 
1.18 
1.47 
1.76 


19.96 
20.77 
21.45 
22.65 
23.82 


0.92 
1.00 


22.51 
23.11 



Table 7. Fiducial points in the plane: (mp6Q6W — m F8l4w). '"F606W- 
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Figure 3. Comparison of three GGCs CMDs (M 92, NGC 6752 and 47 Tuc) with the corresponding isochrones (see text for details). 



and the most up to date in-flight transmission curves of the ACS 
camera. We have named this new in-flight system the WFC/ACS 
Vega-mag photometric system. The proposed calibration is accu- 
rate at the level of a few hundredths of a magnitude. We give the 
basic recipes to calibrate both the observed data and the models, 



and present a set of models by Pietrinferni et al. (2004) already 
transformed into the new ACS Vega-mag system. 

Since the primary purpose of this investigation is to make 
available to the public a complete set of stellar models and 
isochrones transformed to the ACS Vega-mag system, we have de- 
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veloped the following web site where all models and isochrones 
can be retrieved: http:// www.te.astro.it /BASTI /index. php. 

We show CMDs for 5 Galactic globular clusters, spanning a 
metallicity range -2.2 <[Fe/H]< -0.04; these CMDs have been 
obtained from ACS observations, calibrated to the ACS Vega-mag 
system, and used to test our models. We also give the fiducial points 
representing their RGB, SGB, and Main Sequence from the red 
giant branch tip to several magnitudes below the TO. These CMDs 
and fiducial points will be useful for the study of the old stellar 
populations in other stellar systems (Brown et al. 2003). 
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